Chirality and Correlations in Graphene 
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Graphene is described at low-energy by a massless Dirac equation whose eigenstates have definite 
chirality. We show that the tendency of Coulomb interactions in lightly doped graphene to favor 
states with larger net chirality leads to suppressed spin and charge susceptibilities. Our conclusions 
are based on an evaluation of graphene's exchange and random-phase-approximation (RPA) correla- 
tion energies. The suppression is a consequence of the quasiparticle chirality switch which enhances 
quasiparticle velocities near the Dirac point. 
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Introduction — Graphene can be described at low- 
energies by a massless Dirac-fermion (MDF) model [l[ 
with chiral quasiparticles that are responsible for a num- 
ber of unusual properties [2|. The recent experimen- 
tal realization of single-layer graphene sheets 0] has 
made it possible to confirm some theoretical expecta- 
tions, notably an unusual half-quantized quantum Hall 
effect [3, 0, ■ In this Letter we show that quasiparticle 
chirality in weakly doped graphene layers also leads to 
a peculiar suppression of the charge and spin suscepti- 
bilities. We predict that both quantities are suppressed 
by approximately 15% in current samples and that the 
suppression will be larger if uniform samples with much 
lower densities can be realized. At a qualitative level, 
these effects arise from an interaction energy preference 
for MDF states with larger chiral polarization. Our con- 
clusions are based on an evaluation of the exchange and 
RPA correlation energies of uniform spin-polarized MDF 
systems with Coulomb interactions. 
RPA Theory of Graphene — We study the following MDF 
model Hamiltonian (with K = 1), 
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where r 3 is a Pauli matrix that acts on the two- 
degenerate (K, K 1 ) valleys, k is a two-dimensional (2D) 
vector measured from the K and K 1 points, a 1 and a 2 
are Pauli matrices that act on graphene's pseudospin 
(A, B) degrees of freedom, S is the sample area, h" = 

yV^_ i/Sg is the single-spin density operator, and v q — 
2ne 2 /(eq) is the 2D Fourier transform of the Coulomb in- 



teraction potential e 2 /(er). In Eq. the field operator 

C = $k,A>*k,b>*k>,b>*k',a) is a four-component 
spinor. In this work chirality refers to the eigenvalues 
(s,s' = ±) of <r-k/|k|. 

The model ([TJ requires an ultraviolet cutoff, as we dis- 
cuss below. To evaluate the interaction energy we follow 
a familiar strategy Q by combining a coupling constant 
integration expression for the interaction energy valid for 
uniform continuum models. 
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with a fluctuation-dissipation-theorem (FDT) expres- 
sion [7[ for the static structure factor, 
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where n is the total electron density. This form of the 
FDT theorem takes advantage of the smooth behavior 
of the density-density response function along the imag- 
inary axis Xpp'(cM^)- The RPA approximation for the 
interaction energy then follows from the RPA approxi- 
mation for X' 
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where x (q, iSl) is the non-interacting density-density 
response-function. We have derived the following com- 
pact expression for the contribution for an individual 
MDF model channel 
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In Eq. ([5]) ep = vkp where fcp is the channel Fermi momentum, x^ m Eq. ([4]) is constructed by summing the channel 
response function (x MDF ) over valley and spin with appropriate ep values. For a spin- and valley-unpolarized system 
t^(o) _ ^ MDF [with fcp — > fcp = (47rn/g) 1 / 2 ] where g = g s g v = 4 accounts for spin and valley degeneracy. 
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FIG. 1: (Color online) Cut-off A and coupling constant / 
dependence of the regularized exchange energy <5e x in units of 
the Fermi energy ep. The black solid line corresponds to the 
highest value of the cut-off that we have chosen, A = 2.7 x 10 . 



The energy constructed by combining Eqs. ©-© is 
clearly divergent since x increases with q at large q 
and falls only like f2 _1 at large f2. The divergence is ex- 
pected since the energy calculated in this way includes 
the interaction energy of the model's infinite sea of neg- 
ative energy particles. The MDF model can be expected 
to describe only changes in energy with density and spin- 
density at small ep values. For definiteness we choose the 
total energy of undoped graphene (ep = in all channels) 
as our zero of energy. For pedagogical and numerical rea- 
sons it is also helpful to separate the contribution that 
is first order in e 2 , the exchange energy, from the higher 
order contributions conventionally referred to in electron 
gas theory as the correlation energy. Using Eqs. ©-© 
we find that for unpolarized doped graphene the excess 
exchange energy per excess electron is 
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With this regularization the f2 integrals are finite and the 
q integrals have logarithmic ultraviolet divergences. The 
remaining divergences are physical and follow from the 
interaction between electrons near the Fermi energy and 
electrons very far from the Fermi energy as we discuss 
at length later. The best we can do in using the MDF 
model to make predictions relevant to graphene sheets is 
to introduce an ultraviolet cutoff for the wavevector in- 
tegrals, fc c . fc c should be assigned a value corresponding 
to the wavevector range over which the MDF model de- 
scribes graphene. Based on this criterion [9(] we estimate 
that fc c ~ 1/a where a ~ 0.246 nm is graphene's lattice 
constant. The MDF is useful when fc c is much larger than 
fcF in all channels. 



With this regularization the properties of graphene's 
MDF model depend on the dimensionless coupling con- 
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and on A = fc c /fcp. In Eq. {SJ, c ~ 300v is the speed 
of light, a w 1/137 is the fine structure constant, and 
e depends on the dielectric environment of the graphene 
layer. In typical circumstances / ~ 2. A is ~ 10 in the 
most heavily doped samples studied experimentally and 
can in principle be arbitrarily large in lightly doped sys- 
tems. We expect, however, that many of the electronic 
properties of graphene layers will be dominated by dis- 
order when the doping is extremely light. 

In Fig. Q] and Fig. [2] we plot the exchange and correla- 
tion energies of the graphene MDF model as a function of 
/ for a range of A values. Note that both Se x and Se^ FA 
have the same density dependence as ep oc n 1 / 2 apart 
from the weak dependence on A. The exchange energy is 
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positive because our regularization procedure implicitly 
selects the chemical potential of undoped graphene as the 
zero of energy; doping either occupies quasiparticle states 
with positive energies or empties quasiparticlcs with neg- 
ative energies. Note that including the RPA correlation 
energy weakens the A dependence so that the exchange 
energy per electron scales more accurately with ep. It 
is possible to analytically extract the asymptotic behav- 
ior of the exchange and correlation energies at large A 
by Laurent expanding the integrands of Eqs. I®-® in q 
and retaining only the 1/q terms: 

fc x = — fe-p In (A) + regular terms , (9) 
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[Note that £(/ = 0) = 1/3 so that the exchange and 
correlation energies are comparable in size for typical / 
values.] 

Charge and Spin Susceptibilities — In an electron gas, the 
physical observables most directly related to the energy 
are the — ► 0, q — ► charge and spin-susceptibilities, 
normally discussed in terms of dimcnsionless ratios be- 
tween non-interacting and interacting system values. The 
charge susceptibility is the inverse of the thermodynamic 
compressibility K of the system up to a factor of n 2 . For 
the MDF model of doped graphene 
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where &tot includes band, exchange, and correlation con- 
tributions. In Eq. (fl3|) C = (n-[ — nj_)/(n-f + nj_), %s 1 mea- 
sures the stiffness of the system against changes in the 
density of electrons with spin | and spin |. In a 2D elec- 
tron systems the compressibility can be measured [l(| 
capacitively. We note that this type of measurement 
is less difficult when d/i/dn is large as it is in weakly- 
doped graphene. In bulk electronic systems, the spin- 
susceptibility can usually be extracted successfully from 
total magnetic susceptibility measurements, but these 
are likely to be challenging in the case of single-layer 
graphene. In 2D electron systems, however, information 
about the spin-susceptibility can often ll| be extracted 
from weak-field magnetotransport experiments using a 
tilted magnetic field to distinguish spin and orbital re- 
sponse. 

Our results for the charge and spin-susceptibilities are 
summarized in Fig. [3] and Fig. [5] For experiments per- 
formed over the density range over which properties ap- 
pear to be intrinsic in current samples (A between ~ 10 
and ~ 40), these results predict compressibility and sus- 
ceptibility suppression (apparent quasiparticle velocity 
enhancement) by approximately 15%. Both the sign of 
the interaction effect and the similarity of k and \ are 
in remarkable contrast with familiar electron gas behav- 
ior. In 3D and 2D non-relativistic electron gases both 
are @, EH strongly enhanced by interactions, with the 
charge response diverging at intermediate coupling and 
the spin response diverging at very strong coupling. 
Discussion — The qualitative physics of ferromagnetism 
in metals is most transparent at the Hartree-Fock (HF) 
level. Similarly, the mechanism responsible for the un- 
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FIG. 4: (Color online) Cut-off and coupling constant depen- 
dence of the spin susceptibility xs- The color coding is as in 
Figs. [TH 



usual interaction physics of weakly doped graphene be- 
comes clear when the exchange energy is expressed in 
terms of HF theory quasiparticle self-energies. Correla- 
tions do however play an essential quantitative role. For 
doped graphene the contribution of an individual channel 
to the HF theory interaction energy is 
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where s,s' = ± are the chirality indices of the MDF 
bands (i.e. the eigenvalues of the chirality operator de- 
fined above), 
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is the HF self-energy of the undoped MDF model, 
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is the exchange matrix elements between band states 
s',k' and s,k, and 0k, k' is the angle between k and k'. 
For Coulomb interactions, the factor in square brackets 
on the right hand side of Eq. p^|) tends to be larger be- 
tween states in the same band, i.e. states with the same 
chirality. 

The first term on the right-hand side of Eq. (fT4|) is 
similar to the exchange energy of an ordinary 2D electron 
system. Because it is negative and increases with density, 
its contribution to the exchange energy is lowered when 
spins are unequally populated. If this was the only ex- 
change energy contribution, the spin-susceptibility and 



inverse compressibility would be enhanced by interac- 
tions as usual. The unusual behavior comes from the 
second term. For weakly doped graphene it is sufficient 



to expand E k °^ 



to first order in k; £ 
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cally irrelevant constant that is included in the chemical 
potential chosen as the zero of energy by our renormal- 
ization procedure. Expanding to first order in k gives 
the leading interaction contribution to the velocity renor- 
malization [13| of undoped graphene. In agreement with 
previous work we find that for large A and small / 



J_ 

4ff 



MA) 



(17) 



The physical origin of the velocity increase is the loss in 
exchange energy on crossing the Dirac point from states 
that have the same chirality as the occupied negative en- 
ergy sea to states that have the opposite chirality. It is 
easy to verify that this velocity renormalization is respon- 
sible for the leading ln(A) terms in the exchange energy 
and in the exchange contributions to k _1 and Xs ■ ^ ne 
conventional exchange energy contributes negatively to 
k^ 1 and Xg 1 and competes with the Dirac point velocity 
renormalization. 

When correlations are included, the leading ln(A) con- 
tributions to the interaction energy and to kT 1 and Xs 1 ' 
still follow from the (now altered) undoped system quasi- 
particle velocity renormalization. The enhanced quasi- 
particle velocity is tied to the Dirac point, i.e. to the 
switch in chirality, and results in an interaction energy 
that tends to be lower when the chemical potential is 
close to the Dirac point in all channels. Fig. [3] and Fig. [4] 
illustrate RPA theory predictions for experimentally ob- 
servable consequences of the competition between this 
interband effect and conventional intraband correlations. 
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